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Abstract 



Schur process is a time-dependent analog of the Schur measure 
on partitions studied in ^^I- Our first result is that the correlation 
-^ , functions of the Schur process are determinants with a kernel that 

^ I has a nice contour integral representation in terms of the parameters 

1*^ ■ of the process. This general result is then applied to a particular 

f~^ . specialization of the Schur process, namely to random 3-dimensional 

C^ I Young diagrams. The local geometry of a large random 3-dimensional 

diagram is described in terms of a determinantal point process on a 



^^ . 2-dimensional lattice with the incomplete beta function kernel (which 



generalizes the discrete sine kernel) . A brief discussion of the univer- 
sality of this answer concludes the paper. 



a 

> 1 Introduction 

X. 

j^ ■ 1.1 Schur measure and Schur process 

1.1.1 

This paper is a continuation of |16j. The Schur measure, introduced in JE], 
is a measure on partitions A which weights a partitions A proportionally to 
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S\{X) Sx(Y), where s\ is the Schur function and X and Y are two sets of 
variables. 

Here we consider a time-dependent version of the Schur measure which 
we call the Schur process, see Definition |21 This is a measure on sequences 
{A(t)} such that 

Prob({A(t)}) (X nS^*HA(t),A(t+l)), 

where the time-dependent weight S*^*^(A,/i) is a certain generalization of a 
skew Schur function. It is given by a suitably regularized infinite minor of 
a certain Toeplitz matrix. This determinants can be also interpreted as a 
Karlin-McGregor non-intersection probability ^21- The distribution of each 
individual A(t) is then a Schur measure with suitable parameters. 

In this paper, we consider the Schur process only in discrete time, al- 
though a continuous time formulation is also possible, see Section 1^.2. Ill 

The idea to make the Schur measure time dependent was inspired, in part, 
by the paper [TUj which deals with some particular instances of the general 
concept of the Schur processes introduced here. For another development of 
the ideas of ^D], see the paper ^^ which appeared after the results of the 
present paper were obtained. 

1.1.2 

Our main interest in this paper are the correlation functions of the Schur 
process which, by definition, are the probabilities that the random set 

©({A(t)}) = {(t,A(t),-z + i)}, t,zeZ, 

contains a given set f/ C Z x (Z + |). 

Using the same infinite wedge machinery as in ^H] we prove that the 
correlation functions of the Schur process have a determinantal form 

Prob (U C 6({A(t)})) = det (ir(n„ u,))„^_,^,^ , 

and give an explicit contour integral representation for the correlation kernel 
K, see Theorem n This theorem is a generalization of Theorem 2 in |16j . 



1.2 Asymptotics 

1.2.1 

The contour integral representation for correlation kernel is particularly con- 
venient for asymptotic investigations. Only very elementary means, namely 
residue calculus and basic saddle-point analysis are needed to derive the 
asymptotics. 

1.2.2 

We illustrate this in Section El by computing explicitly the correlation func- 
tions asymptotics for 3-dimensional Young diagram in the bulk of their limit 
shape. Three-dimensional diagrams, also known as plane partition, is an old 
subject which has received recently a lot of attention. More specifically, the 
measure on 3D diagrams vr such that 

Prob(7r) oc g'^' , < g < 1 , (1) 

where |7r| is the volume of vr, is a particular specialization of the Schur process 
for which the time parameter has the meaning of an extra spatial dimension. 
As g — > 1, a typical partition, suitably scaled, approaches a limit shape which 
was described in ||3j and can be seen in Figure [7| below. The existence of this 
limit shape and some of its properties were first established by A. Vershik 
|19j . who used direct combinatorial methods. 

It is well known that 3D diagrams are in bijection with certain rhombi 
tiling of the plane, see for example Section 12.3.61 below. Local statistics of 
random domino tihngs have been the subject of intense recent studies, see 
for example IH El UHj and the survey jEj . 

The authors of [H] computed local correlations for domino tilings with 
periodic boundary conditions and conjectured that the same formula holds 
in the thermodynamic limit for domino tilings of more general regions, see 
Conjecture 13.5 in 0. 

From the point of view of statistical mechanics, the argument behind this 
conjecture can be the belief that in the thermodynamic limit and away from 
the boundary the local correlations depend only on macroscopic parameters, 
such as the density of tiles of a given kind. In particular, in the case of 3D 
diagrams the densities of rhombi of each of the 3 possible kinds are uniquely 
fixed by the tilt of the limit shape at the point in question. 



1.2.3 

In this paper, we approach the local geometry of random 3D diagrams using 
the general exact formulas for correlation functions of the Schur process. For 
3D diagrams, these general formulas specialize to contour integrals involving 
the quantum dilogarithm function (j24p . see Corollary ^ 

We compute the q —>■ 1 limit of the correlation kernel explicitly. The result 
turns out to be the discrete incomplete beta kernel, see Theorem |21 One can 
show, see Section 13.1.121 that this kernel is a specialization of the kernel of 
j2] when one of the parameters vanishes in agreement with Conjecture 13.5 
of0. 

The incomplete beta kernel is also a bivariate generalization of the dis- 
crete sine kernel which appeared in jT] in the situation of the Plancherel 
specialization of the Schur measure, see also jTU]. 

1.3 Universality 

Although we focus on one specific asymptotic problem, our methods are both 
very basic and completely general. They should apply, therefore, with little 
or no modification in a much wider variety of situations yielding the same or 
analogous results. In other words, both the methods and the results should 
be universal in a large class of asymptotic problems. 

As explained in Section l3.2.3l this universality should be especially robust 
for the equal time correlations, in which case the discrete sine kernel should 
appear. This is parallel to the situation with random matrices 0. 
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2 Schur process 

2.1 Configurations 

2.1.1 

Recall that a partition is a sequence 

A = (Ai > A2 > A3 > • • ■ > 0) 

of integers such that Aj = for i ^ 0. The zero, or empty, partition is 
denoted by 0. The book ^2] is a most comprehensive reference on partitions 
and symmetric functions. 

Schur process is a measure on sequences 

{A(t)}, teZ, 

where each X(t) is a partition and X{t) = for |t| :§> 0. We call the variable t 
the time, even though it may have a different interpretation in applications. 



2.1.2 

An example of such an object is a plane partition which, by definition, is a 
2- dimensional array of nonnegative numbers 



n 



TT^ 



■i-jj ) 



hJ 



1,2,... 



that are nonincreasing as function of both i and j and such that 



TT 



E 



TCi 



is finite. The plot of the function 

is a 3-dimensional Young diagram with volume |7r|. For example. Figure ^ 
shows the 3D diagram corresponding to the plane partition 

/ 5 3 2 1 \ 

4 3 11 

3 2 1 

V2 1 / 



IT 



(2) 




Figure 1: A 3-diniensional diagram tt 



where the entries that are not shown are zero. 

We associate to tt the sequence {A(t)} of its diagonal shces, that is, the 
sequence of partitions 



A(t) = {ni^t+i) , i > max(0, -t) 



(3) 



It is easy to see that a configuration {A(t)} corresponds to a plane partition 
if and only if it satisfies the conditions 

■ ■ ■ ^ A(-2) -< A(-l) -< A(0) y A(l) y A(2) y ... , (4) 

where X >~ n means that A and fi interlace, that is, 

Al > /Ul > A2 > /i2 > A3 > . . . . 

In particular, the configuration {A(t)} corresponding to the diagram Q is 
(2) ^ (3, 1) ^ (4, 2) ^ (5, 3, 1) y (3, 1) y (2, 1) ^ (1) . 

2.1.3 

The mapping 

Ah^6(A) = {A,-z + |}cZ + | 

is a bijection of the of the set of partitions and the set & of subsets & C Z + | 
such that 

|6\(Z+i)<o| = |(Z+i)<o\6|<oo. 



The mapping 



{A(t)} ^ 6({A(t)}) = {(t, A(t), - ^ + i)} c Z X (Z + i) , 



(5) 



identifies the configurations of the Schur process with certain subsets of Z x 
(Z+ |). In other words, the mapping © makes the Schur process a random 
point field on Z x (Z + i). 

For example, the subset 6({A(t)}) corresponding to the 3D diagram from 
Figure[T]is shown in Figure|21 One can also visualize 6({A(t)}) as a collection 
of nonintersecting paths as in Figure |21 




-5-3-11 3 5 

Figure 2: Point field and nonintersecting paths corresponding to the diagram 
in Figure n 



2.2 Probabilities 

2.2.1 

Schur process, the formal definition of which is given in Definition |21 below, 
is a measure on sequences {A(t)} such that 

Prob({A(t)}) (xnS^*HA(t),A(t + l)), 

where S*^*^(/i, A) is a certain time-dependent transition weight between the 
partition fi and A which will be defined presently. 



The coefficient S*^*^ (/i, A) is a suitably regularized infinite minor of a certain 
Toeplitz matrix. It can be viewed as a generalization of the Jacobi- Trudy 
determinant for a skew Schur function or as a form of Karlin-McGregor non- 
intersection probability [T^ . 

2.2.2 

Let a function 

fcez 

be nonvanishing on the unit circle |2;| = 1 with winding number and geomet- 
ric mean 1. These two conditions mean that log0 is a well defined function 
with mean on the unit circle. For simplicity, we additionally assume that 
(p is analytic in some neighborhood of the unit circle. Many of the results 
below hold under weaker assumptions on as can be seen, for example, by 
an approximation argument. We will not pursue here the greatest analytic 
generality. 

2.2.3 

Given two subsets 

X = {x,}, F = {y,}G6, 

we wish to assign a meaning to the following infinite determinant: 

det(0j,,_xj. (6) 

We will see that, even though there is no canonical way to evaluate this 
determinant, different regularizations differ by a constant which depends 
only on (f) and not on X and Y . 

Since our goal is to define probabilities only up to a constant factor, it is 
clear that different regularizations lead to the same random process. 

2.2.4 

The function admits a Wiener- Hopf factorization of the form 



where the functions 

fceiN 

are analytic and nonvanishing in some neighborhood of the interior (resp., 
exterior) of the unit disk. 

There is a special case when the meaning of (jH)) in unambiguous, namely, 
if (f)~{z) = 1 (or (p^iz) = 1) then the matrix is almost unitriangular and the 
determinant in © is essentially a finite determinant. This determinant is 
then a Jacobi- Trudy determinant for a skew Schur function 



det [(pX-^^+j-i) = sx/M ) ' 

where sx/^{(p'^) is the skew Schur function sx/^ specialized so that 

hk = <l>t, (7) 

where h^ are the complete homogeneous symmetric functions. Note that 

sa/m(0+) = O, /i^A. (8) 

2.2.5 

Definition 1. We define the transition weight by the following formula 

S^{fi,X) = ^s^/^{(f)~)sx/u{(j)'^), (9) 

V 

where sa/,/ (</>"'') is defined by ((7j) and s^/^{(l)~) is the skew Schur function s^/i, 
specialized so that 

hk = 4>Zk ■ 

Note that if the determinant det{(f)Xi~fij+j~i) were unambiguously defined 
and satisfied the Cauchy-Binet formula then it would equal Q because 

(</>-.) = (</'t, )(</>.-.)• 

Also note that because of (jH} the sum in Q is finite, namely, it ranges over 
all z/ such that u G fi, X. 



2.2.6 

In what follows, we will assume that the reader is familiar with the basics 
of the infinite wedge formalism, see for example Chapter 14 of the book JI] 
by V. Kac. An introductory account of this formalism, together with some 
probablistic applications, can be found in the lectures ^Tj. We will use the 
notation conventions of the appendix to |E| summarized for the reader's 
convenience in the appendix to this paper. 
Consider the following vertex operators 



r±(0) = exp ^(log0)=pfc a±k 



vfc=i 



where {log(f))k denotes the coefficient of z'^ in the Laurent expansion of 
\og(f){z). In particular, if the algebra of symmetric functions is specialized as 
in (|7j) then 

(log0)fc = -^. A; = 1,2,3,... , 

where pk is the kth power-sum symmetric function. 

It is well known, see for example Exercise 14.26 in ^I], that the matrix 
coefficient of the operators r± are the skew Schur functions, namely 

{T_{(j))Vf„vx) = sx/f,{(j)^). 

It follows that 

S^{fi,X) = {T.{(j))T4(f))v,,vx). (10) 

It is clear from the commutation relation 

r+(0)r_(0) = e^'=(^°s'^)'=(^°§'^)-'=r_(0)r+(0) (ii) 

that different ordering prescriptions in (jlOj) . which produce different regular- 
izations of (jHl), all differ by a constant independent of fj, and A. 

2.2.7 

Now suppose that for any half-integer m, G Z + ^ we choose, independently, 
a function 

(f)[m]{z) = y^^(f)k[m]z'' 
fcez 

10 



as above so that the series 

^\og (j)[m]{z) 

mall 

converges absolutely and uniformly in some neighborhood of the unit disk. 
This assumption is convenient but can be weakened. The functions (f)[m\ will 
be the parameters of the Schur process. 

Definition 2. The probabilities of the Schur process are given by 
Prob({A(t)}) = ^ n S^N (^(^ - i)' ^(^ + \)) ' 

meZ+l/2 

where the transition weight S^ is defined in Definition Q and Z is the nor- 
malizing factor (partition function) 

Z=^ n S0[m] (A(m- i),A(m+ D) . 

{A(t)}meZ+l/2 

2.2.8 

It follows from (fTUI) that Z is given by the following matrix coefficient 

Z={ W T_{<P[m])T+{<P[m])v^,vA , (12) 

ymeZ+l/2 J 

where W denotes the time-ordered product, that is, the product in which 
operators are ordered from right to left in the increasing time order. 
Using (|TT|) and the following consequence of 



r+(0)t;0 = t;0. (13) 

we compute the matrix coefficient ()12|1 as follows 

Z = expf ^ ^/i;(log0[mi])fc(log0[m2])- 

\m\<m2 k 

Our growth assumptions on the functions (f)[m\ ensure the convergence of Z. 

11 



2.2.9 

It is clear from the vertex-operator description that 

Prob(A(t) = f^) ^ sJl[<P+[m]] sJn 0"[m] j , (14) 

\m<t / \m>t / 

where, for example, the first factor is the image of the Schur function un- 
der the specialization that sets h^ to the coefficient of z'^ in the product of 
(f)'^[m]{z) over m <t. 

This means that the distribution of each individual A(t) is a Schur measure 
[T^ with parameters (fT^ . 

2.2.10 

More generally, the restriction of the Schur process to any subset of times is 
again a Schur process with suitably modified parameters. Specifically, let a 
subset 

{tk}cz, kez, 

be given and consider the restriction of the Schur process to this set, that is, 
consider the process 

{A(A:)} = {A(4)}. 

It follows from the vertex operator description that this is again a Schur 
process with parameters 

#] = n '^N. /,mGZ + i. 

The case of a finite set {tk} is completely analogous and can be dealt with 
formally by allowing infinite values of t^'s. 

2.2.11 

The restriction property from Section 12.2.101 forms a natural basis for con- 
sidering the Schur process in continuous time. 

For this we need a function £(2;, s) of a continuous variable s G M which 
will play the role of the density of log0. The restriction {A(tfc)} of the 



12 



process A(t), t G M, to any discrete set of times {t^} is the Schur process 
with parameters 




m](2;)=exp / L{z,s)ds\ , m e Z 



2 



2.2.12 

Let q G (0, 1) and consider the probabihty measure DJlq on the set of all 3D 
diagrams such that 

Prob(7r) oc gl^l . 

We claim that there exists a particular choice of the parameters of the Schur 

process which yield this measure under the correspondence ^. Concretely, 

set 

f(l-gH^)-i m< 0, ^ , 

hD[m]{z) = \) Vl -1^ -1 ^n rneZ + l (15) 

I (1 — q' 'z ) , m > , 

It is well known that if the algebra of the symmetric functions is specialized 

so that 

/^fc = c^ k = l,2,..., 

for some constant c, then 



"/^ '0, /XT^A. 



Therefore, we have 



S<^3dM (A(m - i), A(m + i)) = 

'^m(|A(m-|)h|A(m+|)|)^ ^ < Q, A(m - i) ^ A(m + |) or 

m > 0, A(m - I) ^ A(m + |) , 
, otherwise . 

It follows that for the specialization (fTK|) the Schur process is supported 
on configurations of the shape (^ and the weight of a configuration (Q is 

proportional to 

qL \Ht)\ = qW . 
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In particular, the partition function Z becomes 

E E — k — 

mi,m2 n=l 

which is the well-known generating function, due to McMahon, for 3D dia- 
grams. 

2.3 Correlation functions 

2.3.1 

Definition 3. Given a subset f/ C Z x (Z + |), define the corresponding 
correlation function by 

p{U) = Prob(t/C 6({A(t)})) . 

These correlation functions depend on the parameters (f)[m\ of the Schur 
process. 

In this section we show that 

p{U) = det{K{u,,u,))^^^^^^^, (16) 

for a certain kernel K which will be computed explicitly. 

2.3.2 

Suppose that 

U = {Ui, . ..,Un} , Ui= (ti, Xj) G Z X (Z + |) , 

and the points Ui are ordered so that 

tl<t2<t3---<tn. 

For convenience, we set to = — C)0. From ()10j) and (jllI) it is clear that 

p{U) = ^{Ruv^,v^), (17) 

14 



where Ru is the following operator 

i?^= nr-(0N)r+(0[m]) n l^-.C n r_(0[m])r+(0H) 



m>t„ i=l..n \ ti-i<m<ti 



2.3.3 

Define the operator 

v^.(t) = Ad I n r+(0[m]) n r-(0H)"' ) ■ V'x , 

\m>t m<t / 

where Ad denotes the action by conjugation, and define the operator ^*(t) 
similarly. Note the ordering of the vertex operators inside the Ad symbol is 
immaterial because the vertex operators commute up to a central element. 
It follows from dTTj), (0, and (HH) that 

Now we apply Wick formula in the following form 

Lemma 1 (Wick formula). Let Ai = J2k'^i,k^k o-nd A* = '^f.a*j^^jjl. 
Then 

(n^.^:^0,^0) = det (KAit^j)) , (18) 

where 

j{AiA*v,j,,v,j,), i>j, 

[-{A*AiV^,v^), Kj. 

Proof. Both sides of (J18j) are linear in Oj^ and a*^, therefore it suffices to 
verify (J18p for some linear basis in the space of possible Aj's and ^^'s. A 
convenient linear basis is formed by the series ()19|) as the parameter z varies. 
Using the canonical anticommutation relation satisfied by ip{z) and ip*{z) 
one then verifies (fTBj) directly. D 



We obtain the formula p6|) with 

K{{ti,Xi),{t2,X2)) = 



{^^,ih)^l^it2)v^,V^) , ti>t2, 



{'^l^{t2)^^,{ti)v$,V^) , ti<t2 
15 



Note that for ti ^ t2 the operators "^xiiti) and "^1^(12) do not, in general, 
anticommute so the time ordering is important. However, for ti = ^2 and 
xi ^ X2, these operators do anticommute, so at equal time the ordering is 
immaterial. 

2.3.4 

A convenient generating function for the kernel K can be obtained as follows. 
Set 

fcez+1/2 fcez+1/2 

Set also 

vl/(t, ^) = Ad I n r+(0[m]) H r_(0[m])-i ) ■ ij{z) , 



\m>t m<t 



and define \&*(t,z) similarly. 
We have from 



Ad(r±(0))-v^(z) =0^(^-^)^(^), 
Ad(r±(0))-^*(z)=0^(^-i)-VW. 

and therefore 

vl/(t, ^) = $(t, z) ij{z) , v[/*(t, ^) = $(t, ^)-i ^(^) , 
where 



Ht,z) = ''->''''")' . (20) 



Finally, it is obvious from the definitions that 



zw 

{ip{z) 4j*{w) VH,, Vii) = — , \z\ > \w\ 

z — w 

\/ zw 

-{ilj*{w)tp{z)v(/i,V(^) = , |2:| < lu'l 

z — w 

Putting it all together, we obtain the following 
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Theorem 1. We have 

where the kernel K is determined by the following generating function 

Kt,^t,{z,w)= Yl ^"^«^""'Al(ti,xi),(t2,a:2)), (21) 

^ y/zW ^{ti,z) 
Z-W ^(t2,w) ' 

Here the function $(t, z) is defined by ()20|). and (j2^ is the expansion of fl22p 
in the region \z\ > \w\ if ti > ^2 (^iT'd \z\ < \w\ if ti < t2- 

2.3.5 

In the special case of the measure DJlq on the 3D diagrams the function (J20|) 
speciahzes to the foUowing function 

$3Z)(t,^) = -F^ rT^^^^^IT' "^^^ + 1- (23) 

llm>max(0,t)V-^ V '^ ) 

Consider the foUowing function 

oo 

{z-,qU = \[{l-q-z). (24) 

n=0 

For various reasons, in particular because of the relation ()28|) . this function 
is sometimes called the quantum dilogarithm function ^. 

It is clear that the function (|23|) has the following expression in terms of 
the quantum dilogarithm 



, i' ^ '-' ; 



^:iD{t,z) 



(gV2-V^; q) 



^ {q^/^z;q). 



-, t<0. 
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2.3.6 

In order to make a better connection with the geometry of 3D diagrams, let 
us introduce a different encoding of diagrams by subsets in the plane. Given 
a plane partition 



7r={7rij), i,j = l,2, 



we set 



6(7r) = {(j - t,7i,, -{i+j- l)/2)} , I, J = 1,2, ... . 



There is a well-known correspondence between 3D diagrams and tilings of 
the plane by rhombi. Namely, the tiles are the images of the faces of the 3D 
diagram under the projection 

(x, y, z) ^ (t, K) = {y-x,z-{x\ y)/2) . (25) 

The tiling corresponding to the diagram in Figure ^ is shown in Figure El 




Figure 3: Horizontal tiles of the tiling corresponding to diagram in Figure Q 



It is clear that under this correspondence the horizontal faces of a 3D 
diagram are mapped to the horizontal tiles and that the positions of the 
horizontal tiles uniquely determine the tiling and the diagram vr. The set 

6(7r) C Zx iZ 

is precisely the set of the centers of the horizontal tiles. It is also clear that 
if {A(t)} corresponds to the diagram vr, then 

(t, h) G 6(7r) ^ (t, X + |t|/2) G e({A(t)}) . 

Theorem n specializes, therefore, to the following statement 



Corollary 1. For any set {(tj,/ij)}, we have 

Prob ({(t„ K)} C 6(7r)) = det [KsdHU, hi), {tj, hj))] 
where the kernel K^d is given by the following formula 

K3D{{tl,h^),{t2,h2)) = 

^ f f 1 ^soiti, z) dzdw 



Hl+1 I,- 1^21-1 



(27ri)2 7|^|=i±, J\^\=l:fe Z'W ^3D{t2, w) ^'»i + ^^-'i2 



(26) 



Here the function ^3£){t,z) is defined by (^Hj). < e ^ 1, and one picks the 
plus sign if ti > ^2 d'^d the negative sign otherwise. 

3 Asymptotics 

3.1 The local shape of a large 3D diagram 

3.1.1 

Our goal in this section is to illustrate how suitable is the formula (j^^ 
for asymptotic analysis. In order to be specific, we work out one concrete 
example, namely the local shape of a 3D diagram distributed according to the 
measure Tig as q ^ 1. Our computations, however, will be of a very abstract 
and general nature and applicable to a much wider variety of specialization. 
The reader will notice that the passage to the asymptotics in (j^^ is so 
straightforward that even the saddle-point analysis is needed in only a very 
weak form, namely as the statement that 

f e^'^(") dx^O, M ^ +00 , (27) 

provided the function S{x) is smooth and ^S{x) < for all but finitely many 
points X G 7. 

3.1.2 

Let q = e"'' and r -^ +0. We begin with the following 
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Lemma 2. We have the following convergence in probability 

rVl^2C(3), r^+0, 
where \tt\ is the volume of a 3D diagram vr sampled from the measure dJlq. 
Proof. First consider the expectation of |7r| 

Similarly, the variance of |7r| behaves like 

d __a 

Var IttI = q—- EIttI = o(r ) , 
dq 

whence Var(r^|7r|) -^ 0, which concludes the proof. D 

3.1.3 

It follows that as r — > +0 the typical 3D diagram vr, scaled by r in all 
directions, approaches the suitably scaled hmit shape for typical 3D diagrams 
of a large volume described in [35- Below we will also see this limit shape 
appear from our calculations. 

We are interested in the r —>■ +0 limiting local structure of vr in the 
neighborhood of various points in the limit shape. In other words, we are 
interested in the limit of the kernel (I^Bl) as 



rti^T, rhi^x, 

where the variables r and x describe the global position on the limit shape, 
in such a way that the relative distances 

At = ti — t2 , Ah = hi — h2 , 

remain fixed. This limit is easy to obtain by a combination of residue calculus 
with saddle-point argument. Since the measure 9Jlg is obviously symmetric 
with respect to the reflection t \-^ —t, we can without loss of generality 
assume that r > in our computations. 
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3.1.4 

We have the following r — > +0 asymptotics: 



\n(z;q)oc^r~^ / — — dw = -r'^ di\og(l - z) , (28) 

Jo w 



where 

dilog(l - z) = ^ - , 1^1 <1 



n^ 



analytically continued with a cut along (1, +cx3). 
Introduce the following function 

S{z- r, x) = -{rl2 + x) In ^ - dilog(l - l/z) + dilog(l - e'^z) 

and recall that we made the assumption that r > 0. The function S{z\ r, x) 
is analytic in the complex plane with cuts along (0, 1) and (e"^, +oo). 
As r ^ +0, the exponentially large term in the integrand in ()26|) is 



e^vr-{S{z-r,x)-S{w-r,x))\ ■ (29) 

The saddle-point method suggests, therefore, to look at the critical points of 
the function S{z;t,x)- Since 

z-^^S{z; r, x) = -r/2 - x - log(l - l/z){l - e'^z) , 

the critical points of S are the roots of the quadratic polynomial 

{l-l/z){l-e~^z) = e-"/2-^. 
The two roots of this polynomial are complex conjugate if 

|gr/2^g-r/2_g-x| <2, (30) 

which can be expressed equivalently as 

-2 In [2 cosh j) < X < -2 In (2 sinh -) . (31) 

In the case when the roots are complex conjugate they lie on the circle 

As we shall see below, the inequahty (jHUj) describes precisely the possible 
values of (r, x) that correspond to the bulk of the limit shape. 
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3.1.5 

The following elementary properties of the function S{z; r, x) 



Siz;T,x) = S{z;t,x), 
S{z; T, x) + Sie^z; r, x) = -(r/2 + x)r , 

imply that on the circle 7 the real part of S is constant, namely, 

^S{z;r,x) = -{r/2 + x)r/2, ze^. 
On 7 we also have 

z—S{z] r, x) = 3r/2 - x - In |e^ - zp , 267, 
dz 

From this it is clear that when the critical points of S are complex con- 
jugate, they are the points of intersection of two following circles 

{z„ zj = {|^| = e^"} n {|z - e^l = e'^/^-^l'} . (32) 

This is illustrated in Figure ID which also shows the vector field 

V{^S{z;t,x)) = ^^S{z), zej. 
e^ dz 



3.1.6 

Now we are prepared to do the asymptotics in (PB|) . We can deform the 
contours of the integration in (j26j) as follows 



KsdHU, hi), {tj, hj)) = / dz dw ... 

where dots stand for the same integrand as in pHjl . < e ^ 1, and we pick 
the plus sign if ti > t2 and the negative sign otherwise. 

Now we define the contours 7>, 7<, 7+, 7_. This definition will be illus- 
trated by Figure El The contour 7> is the circle \z\ = e^^"^ slightly deformed 
in the direction of the gradient of ^RS", see Figure EJ Similarly, the contour 
7< is the same circle \z\ = e^^"^ slightly pushed in the opposite direction. The 



22 




Figure 4: Gradient of ^S{z) on the circle 7 = {|2;| = e^"} 



7< 



7- 




Figure 5: Contours 7> (dashed), 7< (dotted), and 7-1- 
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contours 7-1- are the arcs of the circle \z\ = e"^/^ between z^ and Zc-, oriented 
toward z^. 

Deforming the contours, and picking the residue aX z = w, we obtain 



where 



with the same integrand in as in (PUJ) and 

"(2) 1 f (gi/2+*2w;g)oo dw 



^3D((ti, 


hi), {tj, 


h,)) = 


/•(I) rm 

J ^J 


r_ 


1 


f dz 


f dw ... 


J 


{2my 



2m J^^ (gV2+tl^t;;g)^^Ah+Ai/2+l ' 

where we choose 7+ if ti > ^2 and 7_ otherwise. As we shall see momentarily, 

/ ^ , r -. +0 , 
while j has a simple limit. 

3.1.7 

It is obvious that 

"(2) 1 r dw 

2m L,} ' ^Ah+ At/2+1 



l-e-^w)^'—r--j--, r-^+0. 



7± 

The change of variables w = e'^w' makes this integral a standard incomplete 
beta function integral: 

(2) ^~riAh+At/2) f ^^ 

^ - wr —rrr^rrTTTT , r ^ +0 . 



27ri ./e-r-^± yjAh+At/2+1 



Now notice that the prefactor e '^(^^+^*/2) ^[\\ cancel out of any determinant 
with this kernel, so it can be ignored. We, therefore, make the following 

Definition 4. Introduce the following incomplete beta function kernel 



B±{k,l;z) = ^ f\l-wfw 
2m Jj 



k„,,-l-l 



dw , 



where the path of the integration crosses (0, 1) for the plus sign and (— cxd, 0) 
for the minus sign. 
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3.1.8 



It is also obvious from our construction that the function ^S{z) reaches its 
maximal value on the contour 7< precisely at the points Zc and Zc- The same 
points are the minima of the function S{z) on the contour 7>. The behavior 
of the integral J is thus determined by the term (j29p and, by the basic 
principle (fTfj) . the limit of J vanishes. 

3.1.9 

We can summarize our discussion as follows. Set z^ = e^'^Zc, in other words, 

{z,, z,} = {\z\ = e-^/2} n{\z-l\= e-^l^-^l''} , 53^, > . (33) 

It is convenient to extend the meaning of z^, to denote the point on the circle 
\z\= e""^/^ which is the closest point to the circle |;i; — 1| = e~'^l'^~^l'^ in case 
when the two circles do not intersect. In other words, we complement the 
definition ()33|) by setting 



-r/2 



g-x/2 < gr/4 _ g-r/4 



-x/2 



> ei^ + e 



-t/A 



(34) 



We have established the following 
Theorem 2. Let U = {(tj, hi)} and suppose that as r —>■ +0 

rti^T>0, rhi-^x, « = 1,2, ..., 
in such a way that the differences 



At,; 



ti-t 



J ' 



Ahi 



'H '^j 7 



remain fixed. Then, as r ^ +0, 

Prob{f/ C 6(7r)} ^ det 



B± Ati,-,A/ii,+ 



At- 



where the point z^, = z^:{t, x) is defined in ()33|) anc? fl34|) and the choice of the 
plus sign corresponds to Atij = ti — tj > 0. 



Remark 1. These formulas can be transformed, see Section [3.1.121 into a 
double integral of form considered in jH], Proposition 8.5 and Conjecture 13.5. 
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Remark 2. Observe that the hmit correlation are trivial in the cases covered 
by (|34j) . In other words, they are nontrivial unless the inequality (j3(Jj) is 
satisfied. This means that the inequality (jHUjl describes the values of (r, x) 
that correspond to the bulk of the limit shape. 

3.1.10 

In particular, denote by p* (r, x) the limit of 1-point correlation function 
Kzoi.it, h), (t, h)) -^ p^{t, x) , rt-^T.rh-^x- 

This is the limiting density of the horizontal tiles at the point (t, x). We 
have the following 

Corollary 2. The limiting density of horizontal tiles is 

p*{^,x) = — , 

where 6^, = arg^^, (see Figure{^, that is, 



T e 



^x 



6^ = arccos I cosh — — | . (35) 

The level sets of the density as functions of r and x are plotted in Figure 
ini More precisely. Figure El shows the curves 

/ttt 

9,{t,x) = —, fc = 0,...,8. 

o 

The knowledge of density p^, is equivalent to the knowledge of the limit shape, 
see Sections 13.1.131 and 13. 1.141 We point out, however, that additional analy- 
sis is needed to prove the convergence to the limit shape in a suitable metric 
as in j3]. 

3.1.11 

Corollary [21 can be generalized as follows: 

Corollary 3. The equal time correlations are given by the discrete sine ker- 
nel, that is, if ti = t2 = ■ ■ ■ then 

sm{9^{hi - hj)) 



Prob{t/ C 6(7r)} -^ det 
as r ^ +0. 
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7!-{hi - hj) 




Figure 6: Level sets of the density of horizontal tiles 



3.1.12 

The incomplete beta kernel B± can be transformed into a double integral of 
the form considered in j3] as follows. 
Using the following standard integral 



1 

27r^ 



z " ^ dz 



'\z\=a 1-/5^ 

we can replace the condition 



/3^ k>0, \a(3\ < 1 
-/5^ k<0, \a(3\ > 1 
, otherwise , 



\w - 1\ ^ e-^/^-^/^ 



in the definition of B-|-(A;, /; z^) by an extra integral. We obtain 

1 



B±{kJ;z^ 



^ k 1^ I 1 



(27fiy J jH=e-^/^ 1-z + zw 

|^|=e-/4+x/2 



dzdw . 



(36) 



where the plus sign corresponds to A; > 0. 
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3.1.13 

Let z{t, x) denote the ^-coordinate of the point on the hmit shape corre- 
sponding to the point (r, x) ■ This function can be obtained by integrating 
the density p^: (r, x) as follows. 

Consider a tiling such as the one in the Figure 01 and the corresponding 
3D diagram, which for the tiling in Figure 01 is shown in Figure ^ It is clear 
that the 2;-coordinate of the face corresponding to a given horizontal tile 
equals the number of holes (that is, positions not occupied by a horizontal 
tile) below it. It follows that 

z{r,x)= / (1 -p*(r, s))c?s. (37) 



Here, of course, the lower limit of integration can be any number between 
— oo and the lower boundary of the limit shape given by the equation ()31|). 
Integrating ()35|) we obtain the following formula 

This integral can be evaluated in terms of the dilogarithm function. From 

()25|) we can now compute the other two coordinates as follows 

T T 

xir,x) = z{t,x)-X-^, y(^,x) = ^(^,x) -X+ 2 ' (39) 

which gives a parametrization of the limit shape. A plot of the limit shape 
is shown in Figure 13 

3.1.14 

To make the connection to the parametrization of the limit shape given in 
[S], let us now substitute for p,, in the integral (pTjl the formula (jSHl)- After 
a simple coordinate change we obtain 

1 /" f^'^ du dv 



47r2 J J^ I _^ gx/2+T/4+m _^ gx/2-T/4+TO " 

Integrating this in x we obtain 



z{t, x) = ^ [[ In I 1 + e^/2+r/4+» ^ gX/2-r/4+^^, \ ^^ ^^ 



27r2 



2tt 
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Figure 7: The limit shape 

This together with ()39|) is equivalent to following parametrization of the limit 
shape found in j2j 

ix,y,z) = ifiA,B,C)-2\nAJiAB,C)-2\nB,fiA,B,C)-2\nC), 

where A,B,C > and 

f.27r 



f{A,B,C) 



27r2 



In L4 + 5e*" + Ce™ dudv. 



This parametrization is manifestly symmetric in x, y, and z. However, it in- 
volves integrals that are more complicated to evaluate than the parametriza- 
tion (PH|) . 

Note that the shape considered in jH] differs from our by a factor of 2 due 
to different scaling conventions. 
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3.2 Universality 

3.2.1 

The reader has surely noticed that all that we really needed for the asymp- 
totics is to be able to deform the contours of integration as in Figure so 
that the points of the intersection of 7> and 7< the minima of the function 
^S{z) on one curve and the maxima — on the other. The intersection points 
are then forced to be the critical points of the function S{z). 

It is clear that this principle is very general and applicable in a potentially 
very large variety of situations, such as, for example, in the case of Plancherel 
measure, see below. In particular, the existence of contours of a certain kind 
is a property which is preserved under small perturbations. 

3.2.2 

One such perturbation would be to consider anisotropic partitions. The 
anisotropy in the t-direction is especially easy to introduce: one just should 
replace g'™' in (J15p by g^(™) for some function V. 

3.2.3 

Observe that the asymptotics of the equal time correlations are especially 
easy to obtain in our approach. This is because 

dz $(ti,2:) _ 
nes^=^ -r-T- r — i , ti — 12 , 

Z — W 'P[t2, W) 

and hence the analog of the integral J becomes simply the integral 

<^^ 1 r dw 



2ni J^^ w^^+i ' 
which leads to the discrete sine kernel in the variable Ax. 



3.2.4 

Let us illustrate these general remarks by briefly discussing the asymptotics 
for the poissonized Plancherel measure PJ. It corresponds to the following 



30 



specialization of the Schur process 

>pianch["^J = ■\ e^^~' , m = I , 

, otherwise , 

where a > is the poissonization parameter. In particular, only the partition 
A(0) is nontrivial. The correlation kernel at t = specializes to 

{2myjJ z-w 

which can be expressed in terms of Bessel functions of the argument 2^/a, 
see PEj. 

3.2.5 

The a — * oo asymptotics of the kernel -ft'pianch is easy to obtain from the clas- 
sical asymptotics of the Bessel functions, see P . It is, however, instructive 
to see how this can be done even more quickly in our framework. Assume 

that 

X y 



'a a/o; 
in such a way that A = x — y remains fixed. The critical points of the action 

•S'pianch = Z- Z~^ - -InZ 

are complex precisely when |^| < 2, in which case they are the points e^*^, 
where 

9 = arccos(^/2) . 

So, the same argument as we employed above immediately yields the follow- 
ing formula from pp. 

sing A 
A Planch (a;, y) -^ — T— . 

ttA 
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3.2.6 

Of course, a finer analysis (which was carried out in pP) is needed to justify 
depoissonization in the asymptotics. In our situation, a similar problem is 
to pass from the g — > 1 asymptotics of the measure DJlq to the asymptotics 
of the uniform measures on partitions of a given volume N a.s N ^ oo. In 
other words, further work is needed to verify the equivalence of ensembles in 
the asymptotics. 

3.2.7 

Similarly, finer analysis is needed to work with the asymptotics at the edges 
of the limit shapes where one expects to see the Airy kernel appear. In the 
edge scaling, the following equivalent version of the formula (|22j) 

K{{t,, X,), ih, x,))= f; [;.-+-^-"] 1^ , (40) 

may be helpful because, by analogy to the situation with the Plancherel 
measure [TJ |H], one could expect to see this sum to become the integral 
representing the Airy kernel 



M{x)k^{y)-M{x)ki{y) 



/■oo 

/ Ai(x + s) Ai{y + s) ds . 
Jo 



x-y Jo 

Further discussion of the Airy-type asymptotics of the integrals of the form 
(El, (EHD can be found in HT]. 

3.2.8 

Recently, the techniques of this paper were used in ^ to prove that, indeed, 
the boundary of a random 3D Young diagram converges to the Airy process. 



A Summary of the infinite wedge formulas 

Let the space V be spanned hj k, k & Z+j. The space A^V is, by definition, 
spanned by vectors 

Vs = S1AS2AS3A ... , 
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where S = {si >S2>---}cZ + |is such a subset that both sets 

5+ = 5\(Z<o-|) , 5_ = (Z<o-|)\5 

are finite. We equip A^V with the inner product in which the basis {vs} is 
orthonormaL In particular, we have the vectors 

^'A = Ai - I A A2 - I A A4 - I A . . . , 

where A is a partition. The vector 

V0 = -| A -| A -| A ... 

is called the vacuum vector. 

The operator ipk is the exterior multiplication by k 

i^k{f) = kAf. 

The operator ipl is the adjoint operator. These operators satisfy the canonical 
anti-commutation relations 

all other anticommutators being equal to 0. We have 
The operators q;„ defined by 



a. 



Y^ i^k-ni^l, n = ±1,±2, 



fcez+l 



satisfy the Heisenberg commutation relations 

see the formula. Clearly, a* = «_„. It is clear from definitions that 

[a„,^(z)] = z"V^(z) , [an,r{w)] = -w^riw) (42) 

and also that 

q;„ f = , n < . (43) 
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